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I. INTRODUCTION 



In quantum mechanics, the state of a system is represented by a density operator on a 
Hilbert space. The state can be pure or mixed. To each density operator one can associate 
several quantum phase space distributions which provide a quantum analogue of the classical 
phase space distribution. Prominent among these are the Wigner function 1 , the Q-function 2 
and the P-function 3 distribution each with its own special features. These three quantum 
phase space distributions belong to a one parameter family of quantum phase space distri- 
butions introduced and investigated by Cahill and Glauber 4 and by Agarwal and Wolf 5 ' 6 . 
Over the years, quantum phase space distributions have not only proved to be useful compu- 
tational tools by enabling one to transcribe operator equations into c-number language, but 
have also led to new concepts such as non clasical states of radiation. Recent developments 
in quantum state reconstruction 7 , have made it possible to measure some of the quantum 
phase space distributions directly 8 . They are no longer auxiliary concepts useful only for 
computational purposes but have acquired a meaning in their own right. Several schemes 
for direct measurements of quantum phase space distributions 9 or the positivized versions 
thereof 10 as well as those in the context of atoms 11 have been proposed and are likely to be 
experimentally realized in the years to come. 

In this work, we examine the structure of the quantum phase space distributions from 
the point of view of thermofield dynamics 12,13 . In the conventional formulation of quantum 
mechanics, pure states and mixed states are treated on an unequal footing. The formalism 
of thermofield dynamics overcomes this drawback by doubling the Hilbert space. In this for- 
malism, density operators describing pure or mixed states are represented by a state vector 
in the doubled Hilbert space. The dynamical equations in both cases also acquire the struc- 
ture of a Schrodinger equation even when dissipation is taken into account. The usefulness 
of this formalism for practical purposes can be seen from the exact algebraic solution 14 of a 
class of master equations 15 describing coupled dissipative nonlinear oscillators. It is shown 
here that, in the framework of thermofield dynamics, quantum phase space distributions 
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can be expressed as overlap between the state corresponding to the density operator and 
"thermal" coherent states. To put it differently, quantum phase space distributions appear 
as coefficients of expansion whem the state vector corresponding to a given density operator 
is expanded in terms of "thermal" coherent states. That this aesthetically satisfying picture 
is also useful is demonstrated by giving an exact and straightforward algebraic treatment for 
the evolution of the quantum phase space distributions for the case of the master equation 
describing a dissipative nonlinear oscillator. 

A brief summary of this work is as follows. In sec. II we briefly review the definition 
and properties of quantum phase space distributions. In sec. Ill we outline the formalism of 
thermofield dynamics and show how the statements in sec I translate into the thermofield 
dynamics notation. Some illustrative examples of computation of quantum phase space 
distributions are presented in sec. IV. Finally, in sec. V we apply the machinery developed 
in the preceding sections to obtain a complete picture of the time evolution of the phase 
space distributions corresponding to the density operators which evolve according to some 
standard master equations. 



In a series of papers Agarwal and Wolf have developed a general formulation of quantum 
phase space distributions. Crucial to their formulation is the notion of a A-operator. Here 
we shall confine ourselves to class of A-operators relevant for the present work. Consider 
the family of operators 



II. QUANTUM PHASE SPACE DISTRIBUTIONS 




(a-\)\p\ 2 -(pa*-f3* a ) 



a < 1 



(1) 



where D(a) = exp(aa 



a* a). These operators are hermitian 



(2) 



and have the following properties 
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Tr[A (a) (a, a*)] = n5 2 (a) , 



(3) 



Tr[A^ a \a,a*)A^- a \(3,(3*)} = TT5 2 (a-(3) , 



(4) 



which follow from the fact that Tr[D'(a)D(f3) = n5 2 (a—f3). With the help of these operators 
one can associate with a density operator p a class of quantum phase space distributions as 
follows 



$<f\a,a*) = -Tr[pA {a \a,a*)} 

7T 



(5) 



Conversely, p can be expressed in terms of these in the following manner 

p = J d 2 a& p a) (a, a*)A^- a) (a, a*) . (6) 

For a = l, 1/2, 0, $W (a, a*) respectively correspond to the P-function P(a,a*), the 
Wigner function W(a, a*) and the Q-function Q(o;,a*). The family of quantum phase 
space distributions defined in this way are related to each other as follows 



a — P 



; b < a , 



(7) 



$ p a \a,a*) = exp 



(a-b) 



d 2 



dada* 



(8) 



It may be noted that the class of A-operators may also be expressed in the following manner 



A ia \a,a*) =D{a)p^ ] D\a) , 



(9) 



where 



= - f d 2 (3D{(3)e^> 
ix J 



(10) 



The operator p can be expressed in a more familiar form as 



p ( o ] = (1 " e-")e 



0\-0ata 



1-a ' 



(11) 
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and is easily seen to have the structure of a "thermal" vacuum. It correspnds to a genuine 
thermal vacuum only for a < 0. The family of quantum phase space distributions can 
therefore be compactly defined as follows 

^;\a,an = -Tr[pD(a)p { a) D^a)] . (12) 

71 

This observation will be crucial to later developments. It may be remarked here that, 
expressed in this way, the quantum phase space distributions corresponding to a density 
operator p are identifiable as the generating function of the number distribution of the 
density operator D j( (a)pD(a). 

& a \a,a*) = ] —^Y j \ n <n\D\a) P D{a)\n> ; X — —a/(l — a) . (13) 

That the family of quantum phase space distributions <&^(a,a*) can be expressed as in 
(13), was, to our knowledge, first noticed by Moya-Cessa and Knight 16 . 



III. THERMOFIELD DYNAMICS 

In TFD one associates, with a density operator p acting on a Hilbert space H, a state 
vector |p Q >,l/2<o;<lin the extended Hilbert space H ® H* so that averages of 
operators with respect to p acquire the appearance of a scalar product. 

< A >= TrAp =< p x - a \A\p a > . (14) 

The state | p a > is given by 

\p a >=p a \I> , (15) 

where 

p a = p a ®I , (16) 

and 
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\I >=Y, \ N > ^ \ N >=Y, \ N,N > , (17) 

where | N > constitute any complete orthonormal set in H. The state | J > is simply the 
counterpart of the resolution of the identity 

| / >= ]T | N >< N | , (18) 

in terms of a complete orthonormal set | N > in H. In particular, if 

p\N > =p N \ N > , (19) 
\p a > = T,Pn\ N > N > ■ (2°) 

TV 

(It may be noted that for any density operator the states | p a >,l/2 < a < 1 have 
a finite norm in the extended Hilbert space. This is not necessarily so with the states 
| p 1_a >, 1/2 < a < 1, which include the state | / >. These states are to be regarded as 
formal but extremely useful constructs) 

In this work we shall set a — 1. This representation of the density operator as a state 
vector in the extended Hilbert space is the only useful one for discussing dissipative dynamics 
which will be of interest to us later. In this representation, for any hermitian operator A, one 
has 

< A >= Tr(Ap) =<A\p> . (21) 

Further, for | iV >, we choose the number state | n > and introduce creation and annihilation 
operators a\a\ a, and a as follows 

a | n, m >= \pa \ n — 1, m > , a | n, m >= \pm \ n, m — 1 >, (22) 
a) | n,m >= y/n + 1 | n + 1, m > , a) \ n, m >= y/m + 1 | n, m + 1 > . (23) 

The operators a and commute with a and . It is easily seen that the operators a and 
a) respectively simulate the action of a* and a on | n >< m \ from the right. From the 
expression for | I > in terms of the number states 
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|/>=]T|n,n> , (24) 

n 

it follows that 

a\I>=a)\I> ;a^\I>=a\I> , (25) 
and hence for any operator 

A(a\a) = J2® P , q ^ p a q , (26) 
p,q 

one has 

A | / >= A ] | / > , (27) 

where A is obtained from A by making the replacements (tilde conjugation rules) a — > 
a, a* — > a*, a — > a*. 

With the help of the formal constructs introduced above, it becomes possible to represent 
any density operator p as a state vector in the extended Hilbert space. Thus, for instance, 
the density operator for a coherent state p =\ a >< a | can be represented as | p >= 
D(a)D(a*) | 0, >. Similarly for the density operator corresponding to a thermal state one 
has 

p = (1 - e ^)e~^ a -> | p >= (1 - f)e fK+ | 0, > , (28) 

where / = e _/3 and K + = aJcr. The operator K + together with K_ = ad and K 3 = 
{at a + at a + l)/2 satisfies the algebra of su(l, 1). 

[K^,K + }=2K 3 ; [K 3 ,K ± ] = ±K ± , (29) 

with K Q = {ata — at a) as the Casimir operator. Use of the disentangling theorem for 
s«(l,l) 14 

exp( 1+ K + + 73K3 + 1-K.) = exp(T + K + ) exp((2 log ^7^) exp(r_K_) , (30) 

where 
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± 



27 ± sinh r— ( 20 



20 cosh — 73 sinh ' v y 20 cosh — 73 sinh / 

with 



(31) 



2 = (7 3 2 /4) - 7+7- , (32) 
enable us to write the state | p > in (28) as 

I p >= e n(K++K.-2K 3 J 0?0 > _ (33) 

Turning to dynamics, consider, for example, the master equation for a nonlinear oscillator 
d 1 1 

— p = —i[H, p] H — 7(n + l)(2apa} — a)ap — pa) a) H — r yn(2a'pa — aa) p — paaJ) , (34) 

where H = uja) a + x{<^ a) 2 ■ This master equation has been studied by a number of authors 15 
in the context of nonlinear propagation in a Kerr medium. Applying | I > on (34) from the 
right and using (25) this master equation for p goes over to a Schrodinger like equation for 
the state | p > 

d_ 

dt 

where 



p >= -iH \p> , (35) 



— iH = — iuu(a)a — a) a) — ix[{<^a) 2 — (a)a) 2 ] + -7(71 + l)(2aa — a)a — a)a) 

H — 777,(20*0* — aa} — da 1 ) . (36) 

In terms of the operators K + , K_, K 3 and K , the operator —iH can be written as 

- iH = -i(u - x)K + lin + l)K. + 1 nK + - ( 7 (2n + 1) + 2i X K )K 3 + ]p , (37) 

and hence the solution of (35) as 

1 

I p(t) >= exp( 7 oir + -^ft) exp( 7+ ir + + 73K3 + 1-KJ) | p(0) > , (38) 

where 
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7+ = -ynt ■ 7 _ = 7 (n + l)t ; 73 = -(l(2n + 1) + 2ix# )t ; 7o = - x)t ■ (39) 

Using the disentangling theorem (30), (38) can be written as 

I p(t) >= exp( 7o K + i 7 t) exp(r + ^ + ) exp((2 log ^^3) exp(r_K_) | p(0) > , (40) 

The fact that K + , K_ and K 3 have simple actions on | n,m > enable one to solve (35) and 
hence (34) purely algebraically. Detailed expressions for p m>n (t) and the Q-function for an 
arbitrary initial condition may be found in ref.14. 

In the context of the interaction of a single field mode with a nonthermal and phase 
insensitive environment, Agarwal 18 has considered the following master equation 

d 1 1 

—p = K[apa^ + a^pa-(a^a + -)p-p(a^a+-)} , (41) 
In thermofield dynamics notation this translates into 

^ \ P >=k[K + + K^-2K 3 ] \p> , (42) 

so that 

I pit) >= exp(Kt(K + + K^- 2K 3 )) \ p(0) > . (43) 

To conclude this section, we have seen how various density density operators are rep- 
resented in thermofield dynamics. We have also seen how certain master equations could 
be solved purely algebraically using the thermofield dynamics formalism. In the following 
section we shall show how the quantum phase space distributions are to be calculated in 
thermofield dynamics and how the time evolution of the quantum phase space distributions 
corresponding to the the master equations considered here can be exactly determined. 

IV. QUANTUM PHASE SPACE DISTRIBUTIONS IN THERMOFIELD 

DYNAMICS 

From (12) and (21) it follows that 
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where 



^(a, a*) = - < a, a*; a \ p > , 

P 7T 



a, a*; a > = exp [— a[K + + K_ — 2K Z ]} \ a, a* >, 



(44) 



1 



1 - a 



D(a)D(a*)exp 



1 - a 



0,0 > 



(45) 
(46) 



The quantum phase space distributions can thus be written as overlaps between the state 
| p > and "thermal" coherent states. We now consider some examples to illustrate the 
usefulness of this formula for computation of quantum phase space distributions for some 
standard density operators. 
[1] Coherent states 

Consider the density operator p =\ a >< a \ for a coherent state. The corresponding 
state | p > is | a ,«S >. Substituting this in (44) we obtain 



1 



&;\a,a*) = - -< 0,0 | exp 

F (1 — (1)71 

1 



a 



(1 - a)n 
1 

(1 - a)ix 
1 

(1 - a)n 



< 0, | exp 
a 



1 - a 
a 



exp 



exp 



1 - a 



a — ao 



1 - a 

1 , ,2 

a — ao | 



D{— a)D(—a*) \ ao,a^ > 

| ao — a, a*Q — a* > 

< 0,0 \ ao — a,al — a* > 



(47) 



1 - a 

In the limit a — > 1 one obtains the familiar expression for the P-function for the coherent 
state 



Now, since 



§ {1 \a,a*) = - < a, a*; 1 \ p >= 5 2 (a - a ) ■ 

7T 



< a, a*; a \ ao, a^; b >=< a, a*; a + b \ ao, a* Q > , 



(48) 



(49) 



it follows that 



i 



< a, a*; a \ ao, a^; b > 



(l-(a+6)) 



exp 



1 I 12 

1 a — ao 



l-(a+5) 



nS 2 (a — a ) 



if a + b < 1 
if a + b = 1 



(50) 
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This relation will be used later. 

[2] Squeezed thermal coherent state 19 

The density operator for a squeezed thermal coherent state is 

p = (1 - e ~ l3 )D(ao)S(z) exp(-/3<Ja)& \z)D\a ) . (51) 

The corresponding | p > is given by 

\p>=(l-f)D(a )D(a* )S(z)S(z*)exp\fK + ] | 0, > ; / = . (52) 

Substituting this in (44) we obtain 

$<f\a,a*) = (1 ~ /)(1 ~ A) < 0,0 | e A ^D t (a)^ t («^(«o)^(«S)S(z)5( Z *)e^+ | 0,0 > . 



(53) 

By inserting the resolution of the identity in terms of two mode coherent state after exp(AfT_) 
and before exp(fK + ), writing out the action of these operators on the coherent states, and 
carrying out two of the four integrals that occur we obtain 

= ( 1 ~/)( 1 ~ A ) [fl f ^ e -i(l-A 2 )| 7 | 2 +i(l-A)( T a*- 7 *a)-i(l-/ 2 )|5| 2 
9 TT J TT J TT 

< 7 + a - a , A7* + a* - a* \ S(z)S(z*) \ 5, f5* > . (54) 

Finally, on using the known expression for the matrix element of the squeezing operator 
S(z) between coherent states and carrying out the gaussian integrals, we obtain 

a *) = (1 -/)(!- Reefer (1 - /)(! - A) 

where A = —a/ (1 — a) z = re* e and 

X = [(1 - A/) - (A - f)tanh 2 r] \a-a \ 2 

+ /)(1 - A)tan/ir[(a - a ) 2 e" if? + {a* - a* ) 2 e ie ] . (56) 

In the limit of no squeezing, one recovers the familiar expressions for the quantum phase 
space distributions for a thermal coherent state 
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^ a ^ ( i- ( Ia-/)), exp 



a — ao \ 



(57) 



l-a(l-f) 

With this preparation, we are in a position to translate all the relations in sec. II into 
the language of thermofield dynamics 

An arbitrary state | p > can be in terms of | a, a*; a > as follows 

| p >= j d 2 /3^;\(3,(3*) | (3,(3*;l-a> . (58) 

This relation is the analogue of (6). The quantum phase space distribution thus acquire a 
new meaning as coefficients of expansion when | p > is expressed in terms of the "thermal" 
coherent states | a, a*; 1 — a>. 

Taking the overlap of | p > as given in (58) with | a, a*; a — b > and using (49) and (50) 
one obtains the relation (2) between quantum phase space distributions. 

The relation (8) between the quantum phase space distributions comes about in the 
following way. The coherent state | a, a* > has the following structure 

| a, a* >= exp(a(a < — a) — a* (a — a 1 ")) | 0, > . (59) 

The operators (a* — a) and (a — a*)) commute with each other. Now, from (59) it follows 
that 

d d 
(a^ — a) | a, a* >= — \ a, a* > ; (a — e^)) | a, a* >= - \ a,a* > , (60) 

and hence 

d 2 

[K+ + K_ - 2K 3 ] | a, a* >= -(a ] - a) (a - a 1 ")) | a, a* >= - - | a, a* > . (61) 

The relation (8) therefore is a simple consequence of (44) and of the realization of [K + + 
K- — 2X3] as a differential operator while acting on | a, a* >. 
Finally, the representation (58) also leads to the relation 

Tr(Ap) =nj d 2 /^/?, /?*)$^ a) (A (3*). (62) 
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V. TIME EVOLUTION OF QUANTUM PHASE SPACE DISTRIBUTIONS: 



In sec. Ill we discussed some standard master equations and showed how they could 
be transcribed as Schrodinger like equations. In the cases considered we showed that the 
solution of the Schrodinger like equations thus obtained can be written as (38) with appro- 
priate identification of the coefficients that occur there. In this section we wish to examine 
how transcription enables us to obtain a complete picture of the corresponding phase space 
distributions subject to given initial conditions. Two initial conditions are considered 
[1] Evolution from an initial coherent state | p(0) >=| a ,aQ > 

For this initial condition, it follows from (38), (44) and (45) that 

= - <a,a*;a\ e (7o^o + h*) e (r+^ + ) e ((2io g v ^)^) e (r_K-) | ^ > 

= -e^ <a,a* \ e (-a(*r + +^--2JC3) e (r + Jc + ) e ((2io g vri)Jc 3 ) e (r_A:_) aoe -i W t j a *^t > 

TT 

= -eh* < a ,a* \ e (A^ + ) e ((2io g (i-A)^) e (A^) e (r + x + ) e ((2io gv ^)^) e (r-^) | > 



7T 



(63) 



where A = —a/(l — a). Here in the last step we have used the disentangling theorem (30). 
Repeated use of the commutation relations and the disentangling theorem enable us to bring 
the rhs of (63) into the following standard form 

a*) = i < a, a* \ ^W^o^k^^k.) | ^-a* a ^t > ^ (g4) 

where 

p' _ i _ [(i-A)(i-r + )] 
+ [i-Ar_] ' 

r _ , [(i-r_)(i-Ar + )-Ar 3 ] 

-_ [(l-A)v^] 

is — 



[i-Ar+] • 

In this form the rhs of (64) can easily be evaluated. In the case of the nonlinear oscillator 
r± and T 3 are functions of K and (64) can be written as 
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= - ^ 7 (a*aoe-^) m (aaSe^) n (v /r 3) m+n+1 e {r '+- 1)|Q| "e( r '-- 1 )l Q °l' . (66) 



1 ^ 1 



where it is understood that the K Q in the expressions for Y± and T 3 is replaced by (m — n) . 

When x = 0), i.e. in the case of a linear oscillator, T± and T 3 are constants and (66) 
simplifies to 



$W( a , «*) = ^ e |7t e n-l)|a| a e (T_ -IJIaol^^a'aoe-^+aoS^') . ( 6 7) 



[2] Evolution from a given initial quantum phase space distribution 
An arbitrary initial state can be expanded as 

| p(0) >= J rf 2 a $p(J)(«o,«S) \ a , a* ; 1 - a > . (68) 

For the time evolution of the quantum phase space distributions one therefore obtains 

$%{a,a*) = J d 2 a K( a \a,a*,t ; a , a* , a*) , (69) 

where 

K {a) (a,a*,t ; a ,a* ,0) = -<a,a*;a\ e (-yoK +^t) e (j + K + + 13 K 3+1 .K.) | a ^ a *. (i _ a ) > . 

7T 

(70) 

Using the same manipulations as above one finds that 



e 2 



K^(a,a*,t ; a ,a*,0) = — < a, a* | e ( r "^) e « 21og V r s)^) e (r^-) | ao6 — >* e ^ >, 

7T 



(71) 



where 



n = i-(i-A) 



[(i-r + )(A-r_)-r s ] 



[(A-r_)(i-Ar + )-Ar 3 ]' 
~[(i-r_)(i- Ar+) - Ar 3 i , x 

r - = 1 + (1 - A) (A-r-)i-Ar + )-Arj - (72) 



'1 - A) 2 ^ 



3 



[(A-r_)(i-Ar + )-Ar 3 
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The rhs of (71) can be cast into the form (67) for the case the linear oscillator and in 
the form (66) in the case of the nonlinear oscillator. 

Finally, for the master equation (41), owing to its structure, one has the following inter- 
esting result 

$%(a,an = $% Kt \a,a*) , (73) 

which says that $tyja,a*) sweeps through the entire family of the quantum phase space 
distributions $^ (a, a*) ; b < a associated with the initial density operator 



VI. CONCLUSION 

The observation that the family of quantum phase space distributions considered by 
Agarwal and Wolf have the structure as in (12) and the recognition that this form, in 
thermofield dynamics, translates into an overlap as in (44) constitute the key results of this 
work. These results, in turn, enable us not only to understand various relations between 
quantum phase space distribution in an aesthetically statifying manner, but also prove to be 
extremely useful for computational purposes. The latter aspect of the formalism developed 
here is demonstrated by using it to obtain a full time dependent solution for a family of 
quantum phase space distributions for the master equation describing a nonlinear dissipative 
oscillator. We hope that the results presented here will be useful in the context of wave packet 
dynamics in nonlinear systems in the presence of dissipation, a subject, which has attracted 
considerable attention in recent years. 
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